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Abstract 

We  study  the  scattering  of  plane  waves  by  a  confined  region  of 
dielectric  noise,  using  two  different  approaches.  One  approach  represents 
the  scattering  dielectric  noise  as  a  finite  sample  of  a  homogeneous  random 
process.  The  other  approach  represents  the  dielectric  noise  as  a  locally 
homogeneous  random  process. 
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1,  Introduction 

The  scattering  of  acoustic  and  electromagnetic  radiation  by  random  refi'active 
index  fluctuations  of  the  scattering  mediian  (e.g.,  the  troposphere,  the  ocean) 
is  a  problem  of  considerable  contemporary  interest.  For  brevity,  we  shall  re- 
fer to  these  refractive  index  fluctuations  as  dielectric  noise.  The  purpose 
of  this  paper  is  to  compare  two  approaches  to  the  problem  of  the  scattering 
of  incident  plane  waves  by  a  confined  region  of  dielectric  noise.  In  the  first 
approach,  which  is  the  conventional  one,  the  scattering  dielectric  noise  is  re- 
presented as  a  'section*  of  a  homogeneous  (i.e.,  spatially  stationary)  random 
process.  In  the  second  approach,  the  scattering  dielectric  noise  is  represented 
as  a  locally  homogeneous  random  process.  It  is  thought  that  the  second  approach, 
which  is  based  on  a  recent  generalization  of  the  notion  of  stationarity'- J,  may 
be  found  useful  in  the  applications. 

2.  Scattering  by  a  finite  volume  of  homogeneoiis  dielectric  noise 

Let  the  refractive  index  at  the  point  r  be  denoted  by  the  random  func- 
tion (process)  n(r),  and  let  the  deviation  of  n(r)  from  its  mean  be  denoted  by 

An(?)  -  n(?)  -   <n(?)>    , 

where  the  angular  brackets  indicate  the  operation  of  ensemble  averaging.  By 
definition,  the  dielectric  noise  An(r)  has  mean  zero.  We  assume  that  ^n(r)^^l, 
and  that  Zin(r)  is  second  order  homogeneous,  which  means  that 

(1)         <An(?)An(?')>  -  <(An)2>  C(?-r')  , 

where  the  quantity  ^(/^)   ^  is  the  average  dielectric  noise  power,  and 
C(r-r*)  is  the  correlation  function  of  An(r)» 
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Suppose  now  that  a  scalar  plane  wave  f  exp(ik  •  r)  is  incident  upon  a 

3 

cubical  volume  V  »  d  of  dielectric  noise  described  by  the  homogeneous  process 

An(r)«  (The  wave  vector  k  has  magnitude  k  »  (2n/X},  where  X  is  the  wavelength, 
and  points  in  the  direction  of  incidence.)  Suppose  further  that  the  dielectric 
noise  is  weak,  i.e.,  that  there  is   negligible  multiple  scattering  of  the  in- 
cident radiation.  Then,  to  a  good  approximation,  the  scattered  field  observed 
at  a  point  P  in  the  far  field  of  V  is  given  by  the  real  part  of  the  complex 
quantity 

(2)  f(P)  =  -  ^  k2  e^^  f  e^^-"  An(r)  d?  , 

r2i 

except  for  an  arbitrary  phase '--'.  I'or  convenience,  we  choose  the  origin  of  co- 


ordinates 0  at  the  center  of  V.  The  distance  from  0  to  P  is  denoted  by  R,  and, 

2 

by  hypothesis,  obeys  the  far-field  or  Praxmhofer  condition  RX  »  d  .     The  scat- 

tering  wave  vector  K  is  defined  by 

(3)  K  =  k^-k^  , 

.^ 

where  k  has  magnitude  k  «  (2n/X),  and  points  in  the  direction  from  0  to  P.  (See 
s 

Figure  1.)  In  the  case  of  sound  waves,  ■Uie  field  variable  f  is  taken  to  be  the 
pressure.  With  a  slightmodification,  equation  (2)  can  be  used  to  describe  the 
scattering  of  linearly  polarized  electromagnetic  radiation  with  electric  vector 
in  the  direction  a.  One  simply  takes  f (P)  to  be  the  amplitude  of  the  scattered 
electric  field  measured  in  its  direction  of  polarization,  and  replaces  f  in  the 
right  side  of  (2)  by  f    sin  p,  where  p  is  the  angle  between  a  and  k^  LJ»L  J, 

It  is  always  the  case  in  practice  that  the  scattering  volume  V  can  be 
divided  into  a  large  nxiraber  M  of  subvolumes  which  make  substantially  independent 
contributions  to  f(P).   In  particular  K^  {d/ Ji)   ,  where  the  correlation 
distance  Jt    is  some  reasonable  measure  of  the  distance  over  which  the 
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Figure  I.    The   scattering   configuration    (schematic) 
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process  An(r)  is  appreciably  auto-correlated  j   we  assume  that  d  »  J?  ,  It 
follows  by  a  multivariate  Central  Limit  argument  that  f(P)  is  an  approximately 
Gaussian  random  function  of  P,  in  the  sense  tliat  the  distribution  of  the  m  ran- 
dom variables 

f(P^),  ...  ,  f(P^) 

is  m-variate,  complex  Gaussian,  to  a  good  approximation,  provided  that  m  «  M, 
To  the  extent  that  f(P)  is  a  Gaussian  process,  its  complete  statistical  character- 
ization  is  contained  in  its  covariance  <^f(P)f  (P  )/  •  (The  asterisk  denotes  the 
complex  conjugate.)  By  (l)  and  (2),  we  have 

.2  //A_^2< 

/     I.  A^,y\n    f       I        ^  .^  ,  j, 

e  C(r-r')drdr', 


a)       <f(P)f^(p')>   -   °\  ,  ^   k^  e^^^« 


V  yv 


where  R  is  the  distance  from  0  to  P  ,  another  point  in  the  far-field  of  V,  and 

AR  =  R  -  R  •  As  indicated,  the  variables  of  integration  r  and  r'  separately 

ramge  over  V»  The  scattering  wave  vector  K  is  defined  by  (3)  ,  and  the  scattering 

wave  vector  K  by 

(3a)  K  »  k  -  k^  , 

o   s 

where  k     has  magnitude   (2n/X),  and  points  in  the  direction  from  0  to  P  ,     (See 
s 

Figure  1.)     The  integral  in  the  right  side  of  (I4)   can  be  written  either  as 

(5)  [[      ei'^-<?-?').l(">-''0(?-?.)d?d?.     , 

or  as 

If  An(r)  is  known  to  be  a  Gaussian  process,  then  the  Central  Limit  argument  is 
not  needed,  since  a  linear  transform  of  a  Gaussian  process  is  automatically^ 
Gaussian  process.  However,  in  some  applications,  it  is  suspected  that  An(r) , 
which  may  be  due,  for  exan^^le,  to  turbulent  mixing  of  refractive  index  inhomo- 
geneities,  is  non-Gaussian. 


-  li  - 


T  -^   ■»►   -*   ,^7   •^' 

Introdiicing  new  variables  of  integration  p  =  r-r',  ''(=r   in  the  first  case, 

and  p»r-r,''^  =  r  in  the  second  case,  we  find  that  the  integrals  (5)  and 


(5a)  become 

(6) 
and 


w 


(6a)  f  ,  e«  'P  e^*^-"  '"!  C(?)  d?  d^    . 

I 

The  six-diraensional  regions  of  integration  W  and  W  can  be  simply  described  by- 
giving  their  projections  on  the  (p  ,  '??  )  plane,  (Their  projections  on  the 
(p  >  7  )j  (p_j7tP*  ^"^"^  ^P  '''^K^  planes  have  an  identical  appearance.)  Thus, 

the  parallelogram  in  Figure  2a  shows  the  projection  of  W  on  the  (p  , 7-r^  plane, 

I 

and  the  parallelogram  in  Figure  2b  shows  the  corresponding  projection  of  W  . 

To  approximate  the  integrals  (6)  and  (6a),  we  neke  use  of  the  assumption 
d  »  ^  ,  already  used  in  connection  with  the  Central  limit  argument*  Thus,  C(p) 
will  be  substantially  different  from  zero  only  in  the  strips  of  width  x  shown  in 
Figure  2,  It  is  therefore  a  good  approximation  to  write  (6)  and  (6a)  as 


(7) 
and 


°  .«•?  0(?)  d?  f  ei'^-^'>-^  d^   , 
00  _/T 

(7a)  f  e'^'-P  0(?)d?  f  ei<"  )'^  d^   . 

where,  as  indicated,  the  integration  over  p  runs  from  -co  to  oo  in  all  its  com- 
ponents* Since  C(p)  is  a  correlation  function,  it  is  the  Fourier  transform  of 
a  non-negative  function  $  (k)  ,  known  as  the  (fluct\iation)  spectrum  of  An(r), 
i.e.. 


(8)  C(p)  = 


e-^'P  1  (k)  (^   , 


J 


-00 


-  i,a  - 


A  '^x 


Figure  2.    Projections  of  the   regions    W  and    W' 
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and 
(8a) 


l(k)  = 


/" 


(2n)- 


e^^-P  C(^  dp  . 


-00 


Thus,  introducing  the  form  factor 
(9)      G(k) 


1 


;v 


^J,^  ^    sin(k^d/2)  sin(k  d/2)  sin(l^d/2) 


TTdT?" 


^?^ 


itself  a  correlation  function,  we  finally  have  two  approximations  for  the  co- 
variance  of  f(P),  namely 

2iif^<(An)^> 


(10)      <f(P)f*(p')> 

and 

(10a)     <f(P)f*(p')> 


RR 

2nf^<(An)^> 
J 

RR 


k^  e^^^^  1  it)   G(t?') 


l^ei^^«l(K')G(K-K') 


(The  second  of  these  approximations  differs  from  the  first  only  by  having  ][  (K  ) 
in  place  of  $  (K).)  For  practical  purposes,  it  is  enough  to  require  that  (lO) 
and  (10a)  be  good  approximations  only  when  the  magnitude  of  <^f(P)f  (P  )^  is  an 
appreciable  fraction  of  its  value  for  P  «  P  .  It  follows  from  an  examination  of 
the  form  factor  G(K-K  )  that  for  this  to  be  the  case,  K  and  K  must  both  lie 
within  some  cube  of  volume  V  ^  (2n)7V  in  k-space.  Thus,  the  approximations  (10) 
and  (10a)  are  consistent  provided  that 

(11)      J(K)  ->  i  (k')  ,    when  |K-k'|  <  ^    . 

For  any  reasonable  functional  form  of  f  (k),  (11)  is  valid  if  Kd  »  2ji,  a  condi- 
tion always  met  in  practices  In  fact,  the  weaker  condition  K  K.    »  2n  is  usually 
met  in  practice,  and  is  the  justification  for  the  application  of  inertial  range 
turbulence  theories  to  the  problem  of  scattering  by  turbulence-induced  dielectric 
noise. 
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The  similarity  of  the  approximations  (lO)  and  (lOa)  suggests  that  they 
be  combined  in  the  symmetric  form 

(12)  <f(P)f*(p')>  ~   ^"^^<<f")'>  ...^IkARj^  ,(f.|.,  , 

The  significance  of  this  form  of  the  approximation  will  become  apparent  in  the 
next  section.  It  is  easy  to  derive  from  (12)  the  usual  results  pertaining  to 
the  space  diversity  properties  of  the  scattered  radiation  in  the  far-field  ap- 
proximation* In  fact,  (12)  shows  that  the  pattern  of  space  correlations  of  the 
scattered  radiation  may  be  regarded  as  the  ?^aunhofer  diffraction  pattern  of  the 
scattering  volume  V,  and  has  nothing  to  do  with  the  detailed  structure  of  the 
dielectric  noise.  When  P  =  P  ,  (12)  becomes 

5       2nf2v<(An)^>   , 

<|f(P)|2>  ^   __o_X 1.  k^llK)     , 

R 

which  shows  that  J  (K)  is  proportional  to  the  average  scattered  power  received 

at  a  fixed  point  P. 

Another  derivation  of  (10),  (10a),  and  (12)  can  be  given  by  using  the 

harmonic  analysis  of  the  process  An(r  ).  The  idea  is  to  take  the  approximation 

expand  the  process  An(r)  as  a  stochastic  Fourier-Stieltjes  integral,  and  per- 
form the  integration  with  respect  to  r  before  taking  the  covariance  of  f(P)» 
Thus,  we  begin  by  writing 

(13)  An(?)  -   I    e-^^*^  dN(t^   , 


C 


where  the  Fourier-Stieltjes  integral  (13)  is  meant  in  the  mean-square  sense  L-'-'j 
smd  the  integration  extends  over  all  k-space,  as  indicated  by  the  infinite  limits. 
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It  follows  from  the  homogeneity  of  the  process  An(r)  that  the  increments  of 
the  process  N(k)  obey  the  relation 

(II4)  <dN(k)dN*(k')>   -   <(An)^>  $  (k)  6(k-k')  dk  dk*   , 

where  ^  (k)  is  the  fluctuation  spectrum  defined  by  (8)  and  (8a),  and  5(k-k  )  is 
the  Dirac  delta  function.  Using  (I3)  and  (lli),  we  see  that  the  integral  in  the 
right  side  of  (2)  is 

(15)      ^  f  e^^*^An(?)d?  =. 

where  G(K-k)   is  the  form  factor  defined  by  (9).     Equation  (l5)  shows  that  the 
scattered  signal  is  a  weighted  linear  coiiibination,  actually  a  convolution,  of  the 
(orthogonal)  increments  of  the  process  N(k)»     It  follows  from  (2),   (lU) »  and 

(15)  that  the  covariance  of  f(P)  is 

(16)  <f(P)f*(p')>   -   i^Irk^e^^^   f   I   G(K-&G(K'-t')<dN(k)dN*(k')> 

lin  kR  /-ooj-00 

fV<(An)^>   ^^  ^.^^j^  ^  G(K-k)G(K-k)  I  (k)dk  , 


r   - 1 

dN(k)  i 

'    ei(^-"-?d?    = 

G(K-k)dN(k), 

-00 

Jl 

J-oo 

UiSr' 


r: 


Again,  we  require  that  (l6)  be  a  good  approximation  only  when  the  magnitude  of 

■It  I  ' 

<Cf(P)f  (P  )')>  is  an  appreciable  fraction  of  its  value  for  P  =  P,i.e.,  only 

when  the  functions  G(K-k)  and  G(K-k)  overlap  appreciably.  The  overlap  condition 

is  roughly  that  |K-K  |  <  (2n/d).  Since  it  is  easily  established  that  the  form 

factor  obeys  the  composition  rule 


0 


G(K-k)G(K-k)  dk  =  ^^f-    G(K-K  )   , 
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it  follows  that  if  ^  (k)  is  substantially  constant  in  the  range  JK-K  |  <  (2n/d), 
which  means  that  (11)  is  valid,  then  we  can  write 

(17)        G(K-k)G(K-k)  $  (k)dk  ~  -^^  I  (K)G(K-K  )  -  -^^^1  (K  )G(K-K  ) 


J 


Aa^ 


^  i2n)l  Wk^\  ,(tr)  . 
The  derivation  of  (10),  (lOa),  and  (12)  is  completed  by  substituting  (17)  in  (16)« 

3,  Scattering  by  locally  homogeneous  dielectric  noise 

As  an  alternative  approach  to  the  representation  of  a  confined  region 
of  dielectric  noise  as  a  three-dimensional  sample  of  a  homogeneous  process,  we 
can  incorporate  the  feature  of  being  spatially  confined  into  the  process  An(r) 
itself.  Thus,  guided  by  the  one-dimensional  theory  of  locally  stationary  random 
processes'-  -•,  we  say  that  the  dielectric  noise  is  (second  order)  locally  homo- 
geneous if  the  covariance  of  An(r)  is  of  the  form 

(18)  <An(?)  An(?')>   -  <^|An  (^^V^  C(^?«)  , 

Trfiere<[^(An  (^^^))    J/^  is  a  non-negative  function,  and  C(r-r')  is  a  correla- 
tion function  (i.e«,  a  homogeneous  covariance).     Covariancee  of  the  form  (18)  will 
be  called  locally  homogeneotta  co variances*    For  convenience,  we  introdiKje  the  new 
valuables 

(19)  P*    - 


p      »    r-r' 


The  quantity  <^(An(p  ))  >   has  the  meaning  of  a  sliding  power  factor  which  re- 
normalizes  the  average  instantaneous  power  of  the  process  An(r)  to  a  representa- 
tive local  power  level,  chosen  as  the  power  level  at  the  centroid  of  the  points 
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r  and  r*»  We  assume  that 
(20) 


f    <(An(?^))^  4^        <oo   , 


i.e.,  the  integrated  average  power  of  the  process  Zin(r)  is  finite.  When 
<!l(An(p_j_))  ^  is  taken  to  be  the  positive  constant  <^(An)  ^,  equation  (18) 
reduces  to  the  expression  (l)  for  the  case  of  homogeneous  dielectric  noise. 
In  physical  situations  where  the  average  noise  power  is  a  function  of  position, 
but  is  slowly  varying  with  respect  to  the  correlation  distance  of  the  process, 
it  should  be  an  excellent  approximation  to  write  the  covariance  of  the  process 
in  the  form  (18) •  This  is  typical  of  the  state  of  affairs  when  the  dielectric 
noise  is  produced  by  the  action  of  turbulence.  In  fact,  our  definition  of  local 
homogeneity  stems  from  a  need  to  give  mathematical  form  to  Kolmogoroff 's  qualita- 
tive definition  of  local  homogeneity'--',  wMch  expresses  the  physical  fact  that 
the  turbulent  regime  can  vary  from  one  large  eddy  to  the  next* 

The  question  of  the  existence  of  covariances  of  the  form  (18)  is  non-tri- 
val,  and  is  examined  in  [l].  For  our  purposes,  it  will  suffice  to  give  two  examples. 

a*  The  binormal  covariance 

(21)  exp[^U(pyd)^   exp[I(pyd)^ 

is  locally  homogeneous.  This  follows  at  once  from  the  fact  that 

r 

3-2(r/d)2  3-2(r'/d)^  .  e-U(p+/d)^  ^-ipjd)^ 

2 

and  the  fact  that  exp(-x  )  is  both  non-negative  and  a  correlation  function. 

In  the  case  of  homogeneous  processes  (20)  is  not  valid,  corresponding  to  the 
fact  that  homogeneous  processes  have  the  physically  rnireasonable  property  of 
not  being  spatially  confined* 
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Incidentally,  this  example  shows  that  for  a  process  to  be  locally  homogeneous, 
it  is  not  necessary  for  the  average  noise  power  to  be  slowly  varying  with  re- 
spect to  the  correlation  distance,  since  both  factors  in  (21)  have  approximately 
the  same  width* 

b»  It  follows  from  the  closure  of  the  class  of  covariances  under  multi- 
plication that  the  product  of  the  binormal  covariance  (20)  and  any  correlation 
function  is  locally  homogeneous.  In  particular,  if  we  multiply  the  binoraial  co- 
variance  by  any  Gaussian  coirelation  function  exp[^(p_^/d. )  J,  we  obtain  a  locally 
homogeneous  covariance.  Writing 

-(p./d)^       -(P./di)^     ^         'iPjdf 


e 


we  have 


so  that 


{1/dY  =  (i/d)2  +  (i/d^)2    , 

,'2  2  '2  2 

d^id'^       ,  d^;^d£ 

It  follows  that  the  more  general  form  of  the  binormal  covariance 

(22)        ,MP.M)'    ,-(o.M'f    ,     d-  s  d   . 

is  also  locally  homogeneous. 

Consider  now  the  scattering  of  plane  waves  by  locally  homogeneous  di- 
electric noise  (18),  and  suppose  that  most  of  the  contribution  to  the  integral 
.  ^  <^(An(p^))  ^dp^     comes  from  a  certain  finite  region  V.  (For  example, 
for  the  binormal  covariance  (22),  V  is  a  sphere  of  radius  ->-'d.)  Then  it  is  a 
good  approximation  to  write  the  scattered  field  observed  at  a  point  P  in  the  far 
field  of  V  as 


i- 


(23)         f(P)  -  -  ^  k^  e^^  f   e^-"  An(?)  d?  . 

The  only  difference  between  (23)  and  (2)  is  that  the  integrel  extends  over  all 
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space  instead  of  over  a  finite  scattering  volvune  V,     Correspondingly,  the  co- 
variance  of  f (P)   is  given  by 

<f(P)f*(p")>     .    4?-r  k"*  e^''^     C     C  .^<^-^-^-^\i\n(.t)A.(.?')><^<, 

hnSR  j-coj-oo 

which  is  to  be  compared  with  (Ij),  However,  because  of  the  infinite  limits,  it  is 
now  a  simple  matter  to  evaluate  the  double  integral  exactly.  Substituting  from 
(18) ,  we  find  that  the  double  integral  becomes 

(2U)   I    I   exp[i(^).pJexp[i(tK').?J  <(An(J^))2>  C(pJ  dj^  d?_ 


,CD      _  A  ^,     -,  ,00 


[i(^)-P.]  C(?J<i?.  [  exp[i(tK')-?^<(An(?^))2>dp^, 


exp 
f-oo   ■-       -*  y-co 


where  the  new  variables  p.  and  p  are  defined  by  (19).   By  (8a),  the  first  in- 

o    / K+K  1 
tegral  in  the  right  side  of  (2lj)  is  just  i2ny   |  ( -5—]  •  ^®  '^'^  define  a  new 


form  factor 

H(k) 


<(An(?^))^>d?^ 


'-00 


-1 


.GO      ^ 


f   ei^*PK(An(?^))2>d?,  , 


which  is  the  Fourier  transform  of  a  non-negative  function,  and  hence  a  correla- 
tion function.  The  covariance  of  f(P)  can  then  be  written  as 


2nf2    <(An(?^))^d?^ 
y-00 


<f(p)f*(p')>  -     '-'-   ^, ^,^^^^i{^\  nitt')    , 

which  is  to  be  compared  with  (12).  We  note  that  the  quantity  I  ^<(£Sn(p_^))'>dp^ 

2 

plays  the  role  of  the  quantity  V<(An)  >  in  (12),  as  is  to  be  expected.  The 


-^  -^t 


form  factor  H(K-K  ),  which  describes  the  space  correlations  of  the  received  radia- 
tion, may  be  regarded  as  the  IVaunhofer  diffraction  pattern  associated  with  the 
non-uniform  dielectric  noise  power  <^(An(p^))  )>  .  The  fact  that  a  finite  sample 
of  a  homogeneous  process  can  be  considered  an  approximation  to  a  locally  homo- 
geneous process  is  now  seen  to  be  the  motivation  for  our  choice  of  the  symmetric 
The  quantity  (2U)  is  itself  a  locally  homogeneous  covariance j  see  jl]. 
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combination  (K+K  )/2  in  the  apprcocimation  (12).  Finally,  we  note  that  the 
average  scattered  power  <^|f(P)  |  y   received  at  a  fixed  point  P  is  proportional 
to  J  (K),  as  in  Section  2, 

l4«  Conclusion 

We  have  seen  how  to  express  rigoroiosly  processes  which,  although  not 
homogeneous  in  the  large,  are  locally  homogeneous.  When  plane  waves  are  scat- 
tered by  locally  homogeneous  dielectric  noise  and  observed  in  the  Fraunhofer 
region,  it  is  found  that  the  local  struct\are  of  the  noise  determines  the  average 
scattered  power  received  at  a  fixed  point,  whereas  the  overall  structure  of  the 
noise  determines  the  space  correlations  of  the  radiation  received  at  two  different 
points.  It  would  be  worthwhile  to  try  to  develop  a  theory  for  the  case  where 
the  observation  points  are  not  in  the  Fraunhofer  region,  since,  in  typical  radio 
scattering  experiments,  the  antenna  beams  are  not  narrow  enough  to  meet  the 
Fraunhofer  condition  ^  -'  • 
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